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Guided by the analogy to Mie scattering of light on small particles we show that the propagation 
of a Dirac-electron wave in graphene can be manipulated by a circular gated region acting as a 
quantum dot. Large dots enable electron lensing, while for smaller dots resonant scattering entails 
electron confinement in quasibound states. Forward scattering and Klein tunneling can be almost 
switched off for small dots by a Fano resonance arising from the interference between resonant 
scattering and the background partition. 



I. INTRODUCTION 

A counterintuitive feature of relativistic particles is the 
unimpeded penetration through high and wide potential 
barriers, termed Klein tunneling.^ Rooted in negative en- 
ergy states in the barrier, Klein tunneling is caused in 
graphene by a particularity of the band structured The 
intersection of energy bands at the edge of the Brillouin 
zone leads to a gapless conical energy spectrum in two 
inequivalent valleys. Low-energy quasiparticles are de- 
scribed by a massless Dirac equatiorP and have a pseu- 
dospin which gives the contribution of the two sublat- 
tices of the graphene honeycomb lattice to their make- 
up. Chirality — the projection of the pseudospin on the 
direction of motion — is responsible for the conservation 
of pseudospin and the absence of backscattering for po- 
tentials diagonal in sublattice spaced This allows an ex- 
perimental verification of Klein tunneling in a solid state 
system. Early resistance measurements across a barrier 
indicated Klein tunneling^Kl and in a conclusive experi- 
ment Klein tunneling was revealed in the phase shift of 
the conductance fringes at low magnetic field. 8 

Due to Klein tunneling electrons can usually not be 
confined electrostatically. Thus, designs for circuitry 
taken from traditional semiconductor-based electronics 
cannot be used in graphene and unconventional electronic 
devices are designed on optical analogues This has 
lead to the proposal of a potential step as a lens for prop- 
agating electron beams^ or the experimental implemen- 
tation of the counterpart of an optical fiber cabled 

In the present work we study, guided by the analogy to 
Mie scattering of light on small particles j 12 * 13 ! the scatter- 
ing of a plane electron wave on a circular potential step 
in graphene. The set-up we consider is a plane graphene 
sheet on a gated substrate with a separate circular region 
to tune the potential step by applying a bias. Circu- 
lar dots have mostly been analyzed in etched structures 
with a focus on electron confinement as potential hosts 
for spin qbits.^^ Moreover, single dots and voidd^ as 
well as multiple dots arranged in a corrapS have been 
used to model the scattering by impurities or metallic is- 
lands placed on a graphene sheet. The configuration we 
consider has b een studied for bound states in unbiased 
graphen d 20 * 21 ! and for its ray optical scattering properties 



in single and bilayer graphene! 22 * 23 ! We show that the dot 
operates as an electron switch with the preferred scatter- 
ing angles controlled by its size and applied bias. For 
small dots forward scattering and Klein tunneling can be 
almost completely suppressed due to a Fano resonance^ 
between the background partition and the resonant con- 
tribution to electron scattering. 

II. THEORY 

For a gate potential that is smooth on the scale of the 
lattice constant but sharp on the scale of the de Broglie 
wavelength the low-energy electron dynamics is described 
by the single-valley Dirac-Hamiltonian 

H = -iVa + V6{R-r), (1) 

where R is the radius, V the applied bias of the gated re- 
gion, and a = (a x ,a y ) are Pauli matrices. We use units 
such that H = 1 and the Fermi velocity vp — 1. We 
do not consider effects of disorder or edge roughness of 
the step, which are beyond the present continuum ap- 
proach. Our focus is on the analogy to optics. The dif- 
ference between Dirac's equation, governing a two-spinor 
wavefunction in 2D, and Maxwell's equations in 3D lim- 
its this analogy. In particular Klein tunneling, which is 
also found for a smooth gate potential^ or from a tight- 
binding calculation,^ has no optical equivalent. Known 
phenomena from Mie scattering will appear in graphene 
in new guise to satisfy the absence of backscattering. 
To solve the scattering problem, we expand the incident 
plane wave in eigenfunctions in polar coordinates, 
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(q is the wave vector inside the gated region) term by 
term so that continuity of the wavefunction is satisfied. 
The eigenfunctions in polar coordinates to an energy E — 
ak, where the band index a — 1 denotes the conduction 
and a = — 1 the valence band, read 
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J m Bessel's or Z- 
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Hm^ Hankel's function 



of the first kindP^ The scattering coefficients are given 
by 
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and the transmission coefficients by 

H^ +1 (p)Jm(p) ~ H^\p)J m+1 {p) 



H^ +1 (p)J m (Np) - aa'H^(p)J m+1 (Np) 



(!)/ 



(6) 



(7) 



where we introduced the size parameter p — kR and 
the modulus of the refractive index N = \V — E\/\E\, 
and a' is the band index inside the gated region. These 
coefficients satisfy the relations a_ m = a m -i and b- m = 
aa'b m _i. The electron density is given by n = ip^tp and 
the current by j = ifj^aip where ip = ip m +ijj Ici outside and 
if> = ■0 trans inside the gated region. The far-field radial 
component of the reflected current, which characterizes 
angular scattering, reads 



-rcf 



4 
irkr 



a m' a m[ cos (<y +m + l)(j>) 

+ cos((m — m')<fi)]. (8) 



The scattering efficiency, that is, the scattering cross sec- 
tion divided by the geometric cross section is given by 



Q = - E \< 
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FIG. 1. (color online) Comparison of optical Mie scattering 
and scattering by a gated region in graphene. (a) Scattering 
efficiency Q for Mie scattering for positive refractive index and 
(b) for a gated region in graphene with E < V as a function of 
N and p. In (b) the maxima of the scattering coefficients do, 
ai tt2, and 0,3 are indicated by the black, red, blue and yellow 
lines. Note that Q 7^ in graphene whereas Q = in optical 
Mie scattering for N = 1. The reason is that the refractive 
index in graphene is negative which is already apparent in the 
analogue of Snell's refraction lawP (c) Q for Mie scattering as 
a function of the particle radius a and the inverse wavelength 
A -1 for MgO (complex refractive index N = n + ik below), 
(d) Q in graphene as a function of R and E for V = 1 (N for 
V = 1 below). 



III. RESULTS 

The similar form of the scattering coefficients a m in 
optical Mie scatterin g} 12 ! 13 ! anc j j n graphene should allow 
for the same scattering phenomena. Indeed, the scatter- 
ing efficiency as a function of p and N shows the same 
trends: Q features sharp resonances due to the excitation 
of normal modes for large N and small p and a broad and 
overlapping ripple structure for larger p (cf. Fig. [T£i,b). 

Which of these behaviors is realized depends on the 
actual values of the refractive index. In optics the po- 
larizable excitations of the solid determine N. For in- 
stance for magnesium oxide (Fig. [It), a dielectric with a 
strong phonon resonance, the real (imaginary) part of N 
is large only below (above) the phonon resonance which 



entails ordinary (anomalous) optical resonances! 13 * 28 * 2 ^ 
Over most of the spectrum N is small which leads to a 
broad ripple structure. In stark contrast, iV in graphene 
is tied to the ratio of wavenumbers. A small E leads to 
a large N while E — > V entails N — > 0. Thus, for low 
E very sharp resonances are realized which broaden and 
overlap for larger energies (Fig. Ilji). 

Let us now analyze the scattering by a gated region in 
graphene in more detail. For large radius compared to 
the wavelength of the electron scattering shows features 
known from ray optics. Refraction inside the gated re- 
gion gives rise to two caustics which coalesce in a cusp 
(Fig. [2^,). The caustics have already been studied in Ref. 
1221 including a ray-optical derivation for them. Obvi- 
ously, the boundary of the gated region acts as a lens 
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FIG. 2. (color online) Lensing and particle confinement by the gated region. Density n = tp^ip (top row), current density 
j = ip^aip (middle row) and far field radial component of the reflected current j^ et (bottom row) for (a)-(c) R = 100, E = 0.5, 
V = 1, (d)-(f) R = 2, E = 0.1, V = 1, and (g)-(i) 7? = 4, £ = 0.02813, V = I. (a) For large R and moderate £ refraction 
increases n along two caustics, (b) The electron flow around the dot causes (c) peaked forward scattering. The ao mode has 
(d) maximum electron density in the center (e) two vortices in the current field and, (f) preferred forward scattering. The a\ 
mode has (g) a ring shaped electron density, (h) six vortices in the current field and, (i) three preferred scattering directions. 
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FIG. 3. (color online) Schematic representation of the current 
in the gated region for the ao mode. The angular dependence of 
the current (see Eq. (Ill) leads to j r = on the y-axis and to 
j,/, — on the x-axis. The radial dependence in Eq. ( 11 1 leads 
to a concentric circle where j^, — 0. Around the intersection of 
this circle with the y-axis two vortices are formed which stream 
the electron through the middle of the gated region. 
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FIG. 4. (color online) Depletion of Klein tunneling by Fano resonance, (a) Schematic representation of the vortices in the 
current field: (1) and (3) Vortex pattern of the ao and a\ mode. (2) Enhanced counterrevolving top and bottom vortices of the 
ai mode (red) deplete Klein tunneling. (4) Top and bottom vortices are suppressed and electrons stream through the gated 
region, (b) Far field radial reflected current j r r e for cp = (forward scattering) and tp — ±27r/3 as a function of the energy E 
for a gated region with bias V — 1 and radius R — 4.5. For tp = the lineshape is asymmetric due to interference between the 
ao and ai mode while for (p — ±27r/3 it is Lorentzian. (c) f r ct in the far field as a function of the scattering angle tp and (d) 
current density j = tp' atp in the near field for E = 0.0863, R = 4.5 and V = 1. Marker 2 indicates the suppression of Klein 
tunneling. 



focusing the electron beam. The radiation characteristic, 
given by the radial component of the far field reflected 
current f r ci , shows the absence of backscattering and a 
peaked forward scattering (Fig.[2j:). This is surprising as 
the incident electron impacts non-normally on parts of 
the circular gated region which should give rise to reflec- 
tion. It turns out that in the near field the electron flows 
around the dot (Fig. 2b), which results in the peaked 
forward scattering. 

For small radius and low energy of the incident electron 
scattering resonances appear at specific values of E and 
R where one of the a T m reaches unity. In between the res- 
onances only the ao mode is excited off-resonantly. For 
small p the coefficient a r m gives rise to a series of reso- 
nances for pN — jm.s, where j m ,s is the s th zero of J m . 
For small p (or E) and higher orders m their linewidth 
shrinks (hence, the resonances a m >± are not resolved for 
small E in Fig. fb,d). In the limit E = the resonances 
are located at RV = j m ,s with zero lincwidths, in agree- 
ment wit h the bound states found for a dot in unbiased 
graphenePSEB If only one mode is excited the electron 
density, 

n = Kf Wv) + Jl+ifr)) , (10) 



is radially symmetric (Fig. [2ji,g) . For resonant scattering 
the electron density increases dramatically inside the dot 
(Fig. [2^; note the different intensity scale). This indicates 
temporary particle confinement in the gated region. The 
current field inside the dot is 

J = RJ 2 « (Jm+i(<l r ) + J L{ir)) cos((2m + l)0)e r 

+ | a L| 2 " ( J ?n+i(<l r ) - J m(? r )) sin((2m + l)c£)e . 

(11) 

An analysis of this expression (see Fig. [3] for the ao mode) 
reveals a vortex pattern symmetric to the x-axis which 
dominates the near- field (Fig.|2^,h). The incident wave 
is fed into vortices, which trap the electron. Note that 
the particle is confined in vortices, not by total internal 
reflectionPS The vortex pattern of the mode a m is domi- 
nated by 2(2m + 1) vortices close to the boundary of the 
circle which are reflected by 2m + 1 preferred scattering 
directions in the far field. The radial reflected current 

3 f ~ cos((2m+ IV) + 1 (12) 

shows that for the ao mode only forward scattering is 
favoured (Fig. [2^) , while for higher modes more preferred 
scattering directions emerge (Fig. Si for ax). 
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For small electron energies the ao mode is relatively 
broad compared to the sharp resonances of higher modes. 
Constructive and destructive interference between a reso- 
nant a m mode and the off-resonant ao mode can give rise 
to Fano resonances similar to the effect of the hybridiza- 
tion between a continuum and a discrete level in elec- 
tronic transitions.^ In optics they are routi nely used to 
tailor properties of plasmonic structure d 30 ^ 31 ! and ha ve re- 
cently been identified for Mie scattering by a sphere J22I22I 
where they allow to switch from forward to backward 
scattering. The scattering efficiency reveals no Fano sig- 
natures but the angle dependent scattering which de- 
pends also on the phases of the scattering coefficients en- 
compasses interference effects. Remarkably we find that 
over a Fano resonance a small variation of parameters can 
lead to dramatic changes in the near field, such as an in- 
version of the vortex pattern which causes in the far field 
a change in the preferred scattering direction and even 
leads to a suppression of Klein tunneling and electron 
depletion behind the quantum dot. Figure [4] gives the 
example of the a\ resonance (preferred scattering angles 
tp = 0, ±27r/3), over an off-resonant ao background (pre- 
ferred scattering angle tp = 0). For p = the contribu- 
tions of the ao and a\ mode can have the same amplitude 
and interference causes an asymmetric lineshape of j T r ci 
across the resonance with enhancement above and sup- 
pression below. This can be understood from the vortex 
structure in the near field. Above the resonance destruc- 
tive interference suppresses counterrevolving vortices and 



electrons stream through the center of the gated region. 
Below the resonance the couterrevolving vortices are en- 
hanced and suppress forward scattering (Fig.liLb) . Only 
side scattering remains in the far field (Fig. l4p) and, most 
notably, in the near field the current field shows a deple- 
tion of Klein tunneling at the dot (Fig. |2Ji). 

IV. CONCLUSIONS 

To sum up, we demonstrated that a circular gated re- 
gion in graphene may act as a lens focusing a Dirac elec- 
tron beam. Scattering resonances enable particle confine- 
ment and interference effects may switch forward scatter- 
ing on and off. Our results, which should be verified 
by conduction measurements or by mapping the elec- 
tron density at the gated region exposed to a plane elec- 
tron wave by a scanning-tunneling microscope, suggest 
a strong analogy to light scattering by a sphere. The 
gated region may act as a switch in graphene, which by 
tuning its potential manipulates the preferred scattering 
direction. This opens the possibility for the design of 
graphene-based circuitry by spatially structured electric 
biasing. 
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